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Abstract 

We consider the flavour-changing operators associated with the strong axial and 
trace anomalies. Their short-distance generation through penguin-like diagrams is 
obtained within the QCD external field formalism. Standard-model operator evolution 
exhibits a suppression of anomalous effects in K and B hadronic weak decays. A 
genuine set of dimension-eight AS — 1 operators is also displayed. 
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1 Introduction 



Broken symmetries have often proved as useful as exact symmetries to describe the physical 
world. In this respect, particular attention has been paid to the anomalies of Quantum 
Field Theory Their presence in global axial 2 and scale [3] symmetries of massless 
QED and QCD is notably required by the physics of light (pseudo-)scalar particle decays 
into two photons [I] and hadronic transitions between quarkonium levels 

In hadronic weak decays, the QCD axial and trace anomalies have already been advo- 
cated to understand the unexpectedly large branching fractions for B — ► r]'X s [S] and the 
AI =1/2 selection rule in K — ► tttt [Z][Ej) respectively. However, to our knowledge, no real 
attempt has been made so far to generate both pseudoscalar and scalar anomalous glu- 
onic configurations from short-distance QCD corrections to a four-quark flavour-changing 
effective interaction. 

At lowest order in the momentum expansion, the local anomalous operators in request 
consist of a flavour-changing quark bilinear coupled to the colour-singlet gluon density 
G^Gnv in the axial anomaly case and G^G^ in the trace anomaly one. At the one- 
loop level, they are thus generated by gluonic penguin-like diagrams with a heavy virtual 
flavour. In this paper, we infer the corresponding QCD evolution of an arbitrary four- 
quark operator from the propagator of the heavy intermediate flavour plunged in external 
gluon fields and expanded in inverse powers of the mass. The contributions of the leading 
anomalous operators are then displayed. 

Due to our assumptions of a four-quark initial operator and a heavy intermediate 
flavour, the method turns out to be particularly well-suited to the estimation of short- 
distance charm effects in Kaon decays within the Standard Model. Yet, the insight gained 
allows us to clarify the role of heavy quark generated anomalies in B decays too. Besides, 
the formalism developed here is applicable to any four-quark effective interaction induced 
by physics beyond the Standard Model. 

A by-product of our method is the consistent derivation of a complete set of dimension- 
eight operators in standard-model hadronic K decays. Contrary to ref.jH], the central role 
granted here to the heavy flavour propagator immersed in QCD external fields allows us 
indeed to avoid a misuse of the quark classical equations of motion. 



2 Heavy quark induced penguin-like operators 

Let us start with the computation of the one- loop penguin-like short-distance QCD cor- 
rections to the generic four-quark operator 

Q = (q 1 T A q h )(q h T B q 2 ), (1) 

under the assumption of a heavy intermediate flavour q^. The symbols T A,B stand for 
products of Dirac matrices. The quark bilinears are colour-singlets. 

In the path integral formalism, these corrections may be obtained by functional in- 
tegration of the heavy mode in the presence of a classical gluonic background [TUj : 

expi J d 4 x£ ef f = J Vq h T>q h expi J d 4 x (q h pq h - Mq h q h + g w Q) (2) 
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with g w , the effective coupling associated with the weak operator Q, and M, the heavy 
flavour mass. The covariant derivative acting on spinor fields is defined by 

= *fy + g s A^ (3) 

with = A^X a /2, the external gluon field and g s , the QCD coupling constant. The 
Gell-Mann matrices A a are normalized such that tr\ a X b = 2S ab (a, b = 1, 8). 

The Gaussian integral (J2J) is readily performed. It is given by the functional determi- 
nant of the matrix 

A% = [rJ - M S* + g w F B ql x ^ x r A ) 5^ (x-y). (4) 

Colour indices are denoted by i, j = 1,2,3) and spinor indices are understood. Using 
the identity det A = exp tr In A and keeping only the 0{g w ) term, we obtain the (non-local) 
effective Lagrangian £- 9 ™jf- 



J *x£Q f = -ig w J d'x qi x T A 



(5) 



Consequently, the central object needed to compute the local anomalous operators in 
request is the propagator in external field, (P — M)~y = (x (P — M) _1 y\ at x = y. 
We will need the following formal identity: 



P-M)„ I (2lT ft+f>*-M 



Intuitively, it results from the minimal substitution p^ —* p^ + P^ in the momentum-space 
free-quark propagator. The above expression is in fact properly defined by the infinite 
series: 
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P - M 



with 



Sn= j 7^4^71 i-n^j) • (8) 



(2tt) 4 ^-MV V ^ , i>-M / 

Let us emphasize that the derivatives in S n will only affect the involved gluon fields after 
momentum integration, not the external states gi^. Indeed, the 'plunged' propagator is a 
function, not an operator. Eq.0 will thus eventually turn into a well defined expansion 
in 1/M provided the gluonic background fluctuations are small. 

We will now compute the first few S n 's in the gauge-invariant dimensional regulariza- 
tion scheme. The above integrals can be cast into the form 

d d p fi e N n 

(27r) d ( P 2 -M 2 y 



with e = 4 — d and /j,, the regularization scale. The numerators N n are written in a form 
suitable for integration using exclusively the Dirac matrices anticommutation relations 
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{if^lu} = ^9fA,u- In the cases n = 5 and n = 4, we also make use of the following 
decompositions of 7-strings on the Clifford basis after momentum integration: 

7/ll7/127/43 = 5/11/127/13 ~~ 544i/4 3 7/42 + 5/42/437/41 + * £ /Ltl/12M3/37 75 
7/41 7/42 7/43 7/i4 = 5/41/425/43/14 ~~ 5/41/43 5/42/44 "I 5/11/44 5/42/43 ~~ Z£ /41/42/43A4475 
+ * (5/4l/l 2 f7 /44/i3 + 5/4l/t 3 °7i 2 /44 "I" 5/41/44 C/43/I2 

"^5/42/13 °M4/11 + 5/12/44 + 5a43A44 C7 /42A41 ) (10) 

where 75 = »7 7 7 2 7 3 , <r Ml/ = | [7^, 7^] and e Ml /p CT is the Levi-Civita tensor with e 0123 = +1. 
These conventions agree with ref.|12j. 

Gauge invariance requires £1 to vanish. For the next few terms, we easily obtain: 

1 . . /„ , M s 



—7 1 

C _ p p p p pMl/42/13M4 

2 (4tt) 2 M w w M "* 4 

C _ * p p p p p /■■pMl/i2/13/44/45 1 p/41/12/43/44/45\ 

°5 ~~ 2 (4tt) 2 m 2 1 2 3 ^ 5 5 ' ' 

with iV e = I + In 4-7T — 7, and 

2 2 

p/4lM2/13/44 ^/ll/13^/42/44 _| ^/4lM4^/42/13 _ ^/41/12/13/M^ 

+ ig&^aWIM — ^2/43^/41^4 _|_ ^442/44^/11/43 
_ 1^1^2^-/43/44 _ 1/41/44^-/12/13 _ Iq/13/14^/41/42 

3 3 3 

IS 17 8 

pMl«2M3M4M5 1^/12 /43^/44/15^,/41 _| _^M2/44 g/13/45 yil _gM2/45^/43/14^,/41 

30 30 30 

17 18 17 

_| L^MlW^/44/45^,/42 IgMl/M g/43/15 /12 _| L^/41/45^/13/44^,/42 

30 30 30 

8 2 18 

^1^2^/44/45 ^,/43 _| ^Ml/44^/12/15^,/13 1^/41/45^/12/44^/43 

30 30 30 

17 18 17 

_| L^/ 1 lM2^/43/45^,/44 1^/11/13 ^/42/15^,/14 _| _^/41 /15^/12 443 y*4 

30 30 30 

13 17 8 

lg/n/12g/13/44y*5 _| !_^/ t l/ t 3^M2/14^,/15 _gMlM4g/"2M3 yiS 

30 30 30 

p/41 /42/13/14/45 _ _|_lg/41/*2 £ /43/44/45/3^,y 5 _ l^l/*3 e /42/44/15/3^^ 5 _|_ 1^/41 m, £ H2H3HBPjpy 5 

3 2 6 

- ^ WM5 ^ 2M3M4/3 7/375 + ^5 /12M3 ^ lWAt5/3 7/375 + ^5 M2M5 e WAt3M4/3 7/375 

+ l 5 /13/44 £ /41/42/45/3 7/3 75 _ l g M3M5 £ mM2/M/3 7 ^ 75 + 1^/44/45 g/H/42/43/3^^ _ ( 12 ) 

6 2 3 
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We now introduce the conventional definitions for the gluon field-strength tensor G^ u = 
d^A u — d u A^ — ig s [A^, A v \ and its dual G^ v = \e^ V af3G a ^. After some rearrangements 
based on the identities [P^^Py] = igsG^ and £ ll uafdP a P^ = WsG^y, we come to the 
important result: 

5 2 = (n £ - In ¥p\ G^ 

2 (4vr) 2 V M / 

5 3 = -^fiV £ -ln^) VGyst 

5 4 = -7-^2 ^ f G^G"" + -aG^G^s - 3iG£G- a (v, - ----A^G^') 

6 (47r) M \ A Zg s J 

s " = 0fJP {? |D ° G ""' G "" 1 7 " + il ! |0 "» D " G " 1 7 " 

+j {o°G„<„ G""} 7,75 + | {g^, D'-d""} 7„75 - I ^0 < "D a D"G^7") (13) 

with iD^F = [Pu, F] for any field F transforming like the SU(3)c adjoint representation. 

The dominant penguin-induced corrections to the initial four-quark operator (^Q) are 
now explicitly obtained: 

5 

£9 eff = ^ A S n T B q 2 + O (1/M 3 ) . (14) 

n=2 

Remarkably, trace and axial anomalous operators are in principle already generated at the 
lowest possible order in the momentum expansion, as inferred from the colour-singlet part 
of the scalar and pseudoscalar gluon densities contained in . 

Let us now apply these results to the case of the standard-model short-distance weak 
operator evolution. 



3 Anomalies from standard-model operator evolution 

Tree-level integration of the W gauge boson in V — A charged current interactions leads 
to effective operators of the type with T A ^ = 7^ (1 — 75). 

Let us first focus on the well-suited case of charm quark loop contributions to Kaon 
decays. The effective coupling associated with the tree-level operator defined at the M\y 
scale 

Qf =4(d Ll „c L )(c L ^s L ) (15) 

is given by 

g S w M = -V cs V c * d G F /V2, (16) 

with the notation qi, = ^ (1 — 75) q. 

The chiral structure of this operator obviously prevents the even n contributions from 
appearing. The first non-vanishing terms are thus obtained for n = 3 and n = 5. They 
correspond to dimension-six and -eight operators, respectively: 

= £-6d + ^8d + ••• (17) 
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Still, it is worthwhile to analyse them since, as we shall see, they can also carry some 
strong anomaly effects. 



Dimension-six operators 

The dominant operators in the m' 1 expansion ((17)1 have been computed a long time ago 
|13j . Yet, we find it useful to illustrate how the external field formalism works in a simple 
case. It also allows us to fix our conventions. The effect of incomplete GIM mechanism 
[Tl| above m c is at the next-to-leading level ^T5\, i.e. beyond the scope of this paper. So, 

(c) 

the leading short-distance evolution of Q 2 from M\y down to m c only involves the usual 
current-current operators |16j . which disappear once the charm quark is integrated out. 
Inserting the explicit form of S3 in eq. (|14j) . we obtain then at the leading order in the 
strong coupling constant: 

C 6d = ^=V cs V; d -^ ^ Mr} d L D"G^s L (18) 

with [i < m c . Recall that the standard computation of dimension-six gluonic penguin 
operators involves the use of the classical equations of motion: 

D»G,„ = g s Y, (^Y Q )y- (19) 

q=u,d,s 

Applying these equations on Lagrangian (fl8|) together with the identity (A a )y (A a ) M = 
2 (Su5kj — 5jj<5fcz/3) and the relevant Fierz reorderings, we indeed recover the well-known 
current-current and density-density colour-singlet operators: 



k=3 



with 



and 



Qf = A(d L ^s L )(q L ^q L ) Qf = ±(d L y SL )(q R ^q R ) 
Qf = 4(d LW )(q L ^s L ) Qf = -8(d L q R )(q R s L ) 



.(6) _ (6) _ _ 1 (6) _ (6) _ 1 

— <-c: — „ i <^/l — w; 
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In our conventions, q R = i (1 + 75) q and the light flavours q = u,d,s are summed over. 
Dimension-eight operators 

The dimension-eight penguin-like operators are most easily obtained from the contribution 
of 55 to ea. ((TH) . which gives: 
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with 



of 

Q? 
Qf 
Q? 
Qf 



id L [D^G^GnivSL 
id L [G m ,D a G^] lv s L 

d L [G^D a G^) lu s L 
~d L D a D a D v G^s L 



and 



.(8) 



-■, C, 



(8) 



A r (8) 

15' 3 



1 r (8) 



"3' 5 



"15 5 * 



It may be remarked that this effective Lagrangian is invariant under the CPS symmetry, 
as it should [T7] . 

Dimension-eight operators for Kaon decays have already been considered by the au- 
thors of ref.jHj. Yet, their early use of the equations of motion introduces an explicit 
dependence on the light quark masses m^, which is artificial as the 'plunged' propagator 
is fundamentally a function, not a differential operator acting on external spinors. Besides, 
their ensuing truncated expansion in md, s is not allowed. 

Notice that the last operator in eq. (|21|) is still linear in the gluon field-strength tensor. 
As such, it contributes to the finite part of the standard dsg one-loop diagram, providing 
thus an interesting check of our calculation. Indeed, the gluon momentum expansion of 
the 'vacuum polarization' integral that naturally arises from the four-quark approximation 



dx x(l — x) In 



m: 



x(l — x)q" 



1 



(3 



1 



- Mn-f --^ +0(qym 



5 m 



(22) 



displays the correct relative weight between the operators of ea. (|18JI and ea. ()21j) with 
k = 5. 

Let us now investigate the anomalous content of Lagrangian H'21|) . Operators possibly 
involving Ga"G® u and Ga U G^ u gluon densities are necessarily contained in the combina- 
tions 

(23) 



and 



Qf + Q { ! ] =d L D a \G^ ai G^ 



(24) 



respectively. 

Remarkably, the gluonic part of the first operator is a pure SU(3)c octet, and con- 
sequently cannot carry any trace anomalous effect, in contradistinction to the results of 
ref.0. 

On the other hand, the second operator does contain a SU(3)c singlet part. Indeed, 
dropping a harmless total derivative, ea.()24|) can be rewritten: 



id L P^ lv \ '—, '— \ S L + id L { '—, — \ lv P a S L 



X a X b 



2 ' 2 



( X a X 1 



2 ' 2 



G a m G^ 



(25) 
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with P' a = —id a + g s A a , the derivative acting on the left. Taking then the trace over 
colour and Lorentz indices, which amounts to the substitution 

_ sab zv _ 

G a , a G^ b - °—°-fG c pa G?, (26) 

(S) (8) 

and including the coefficients C3 and c\ , we come to the following axial anomalous 
operator: 

G%G^. (27) 
The quark classical equations of motion 

Ps = m s s, ~df f = m d d (28) 

may eventually be used to put this result in a more conventional form. Yet, the matrix 
element of G^ u Ga U between the vacuum and a neutral pion state requires isospin violation 
in the factorization approximation. The charm-induced axial anomaly effect on hadronic 
K decays is thus negligible as far as the AI = 1/2 amplitude is concerned. 



o {8) 
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d L Ps L -d L Pls L 



Dimension-five and -seven operators 

We have seen that the chiral structure of the tree-level operator (|15|) prevents any S^n 
contribution from surviving at the effective Lagrangian level. However, it is well-known 
that dimension-five 'chromo-magnetic' penguin operators do arise in the Standard Model 
once we go beyond the four-quark approximation. A full short-distance calculation with 
massive W propagation gives indeed for both the charm and top contributions |18j : 

C 5d = ^VjjQ^cfqf® (q' = c, t) (29) 

with 

Q(5) = m s d L G^ v s R + m d d R G^ u s L 

and 

cf (Try « M W ) = O (^-] , cf (m q , » M w ) = O (1) . 

So, the longitudinal part of the W propagator induces a chiral structure compatible with 
the S211 operators once the light quark equations of motion are used. However, the price 
to pay is a suppression factor of the order of m s m c /M^ or m s /mt. This can be seen 
explicitly from eqs.(J2U) and (|T3)) in the case of S2. Still, similar suppressions are expected 
to occur also for higher S^ n structures in hadronic K decays. Indeed, the helicity flip 
undergone by external spinors implies a linear dependence on the light quark masses, 
while the denominators can be inferred from the limit Mw — > 00 in the charm case and 
from dimensional arguments in the top case. 

Let us now turn to hadronic B decays. In exactly the same way, the virtual top quark 
contribution is suppressed by a factor of the order of m^/mt- On the other hand, the 
charm loop contribution induced by the tree-level operator Q2 = 4(sl7 aj cl)(cl7 m 6l) is 
rather troublesome since the momenta associated with gluon fields can now be larger than 
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m c . For illustration, the time-like squared momentum of the single gluon configuration 
generated by the 'vacuum polarization' integral (|22|) is likely in the range of < 
q 2 < W/2 for two-body exclusive B decays [TJ2|- If such is the case, the expansion in 
inverse powers of the charm mass obviously breaks down, non-local operators survive and 
an absorptive component may even arise for q 2 > 4m 2 . Consequently, we cannot exclude 
non-negligible anomalous effects associated with S2 n +i chiral structures in B — > r/'K, 
unless the factorization approximation is called for help. 



Factorized amplitudes 

The factorized part of the B — ► rj ' K decay amplitude reads: 



rfK Q 



(c) 



B 



-3 W |c7^75c| 0) (K\sj%\B) 



(30) 



The fact that heavy quarks can only contribute to light pseudoscalar decay constants when 
propagating in a loop implies: 



<V|c7 M 75c| 0} = (rf 



-iTr 



(31) 



with the trace taken over both spinor and colour indices. Inserting the first few terms of 
the charm quark propagator expansion (|13|) in the above matrix element, we come to the 
conclusion that only part of S5 survives to give: 



(7/ |e/75c| 0) = ( rf 



as 



1 



2 (D a G l3a ) a G%> + G% a (D a G^ 



0) + O(l/O. (32) 



127T m 2 

While this quantity remains to be estimated, it is already interesting to notice that the 
hadronization of the cc pair into the rj' meson proceeds indeed through S2n+i gluon struc- 
tures, unlike what emerges from ref. |20| where the charm quark equation of motion is used 
from the start. 

The axial anomaly part of ea. (|32|) is given by: 



{rf |c7 M 7 5 c| 0) 



1 



AA 



with 



Je)AA 
Jr]> 



32m 2 c 



32m 2 c 



V 



7T V 



. t (c)AA u 



rf 



V 



7T 







-2 MeV. 



(33) 



(34) 



The QCD sum rules calculation of ref.|25 ^ as been used to estimate the last matrix 
element. Consequently, the charm quark induced axial anomaly can only account for a few 
per cent of the experimental B — > rf K decay amplitude in the factorization approximation. 

Note that charm quark loop effects in factorized B — > /o(cr)if decay amplitudes may 
be treated in a similar way. Yet these appear to be even more suppressed: 



(f (a) \c^c\ 0) = (/ (a) 



-iTr 



CP - m o)xx V 1 0)=O(l/mi). 



(35) 
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4 Anomalies beyond standard-model operator evolution 

The minimal non-standard operator leading to S^n-related anomalous effects reads 

Q = (AlCr)(c l s r ) (36) 

in the case of hadronic K decays. This tree-level operator generically arises in multi- 
Higgs models |22| . In this framework, the following dimension-seven operators are directly 
obtained from S4: 

r - n as 1 W^n^ ci7\ 
c 7d - -9wi7T-— 2^ C k Qk 1 37 ) 



with 



and 



2 Air m r 

c k=i 



Q? = d L G^G^s R 

QP = id L G^ a a^s R 
Q { P = d L D a D a G^s R 



JX) _ ! If) _ 3 If) _ _o 17) _ -1 1 
-1 — ±1 c 2 ~ 2' 3 ~~ °> c 4 ~ 2 « • 



In particular, their anomalous content is given by: 

Qaa+i* = \dLSR (g%G^ + ^G^Cr) (38) 

(7) (7) 

if the coefficients c\ and are included. However, the effective coupling constant g w is 
in principle suppressed by the characteristic scale of the new physics involved. 

In the factorization approximation, the trace and axial anomalous operators merged 
in ea.(|38|) are the only contributions of ea.(|37|) to hadronic K decay amplitudes. They can 
thus be directly inferred from the factorization of Q matrix elements. Indeed, the charm 
quark loop contribution in 

tt+tT \Q\ K 0X ) f = -L (tt + it- \cc\ o) (o \d L s R \ K°) (39) 

(tTtt \Q\ K-) f = -L (tt \c l5 c\ 0) <tt- \d LSR \ K-) (40) 

is readily estimated using the 'plunged' propagator, provided we make the replacement 
y — > 1 in ea.l|35|l and 7^75 — ► 75 in ea. (|3*T]) . Let us point out that this simply amounts 
to the well-known substitutions E 



m c cc - -IZLG^G? + 0{f /ml) (41) 
m c ci l5 c - gG^Gf + 0(^/ml) (42) 

for heavy quark contributions to light hadron matrix elements, providing thus another 
check of our central result (|13JI . 
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5 Conclusion 



We have proposed a simple way to systematically generate penguin-like operators from the 
inverse mass expansion of the heavy quark propagator plunged in QCD external fields. 
This expansion has been explicitly performed up to the fifth order. The result, expressed 
in a compact and comprehensible form, is also suitable for other applications like, as we 
have seen, the heavy quark contribution to light meson decay constants. 

This formalism has allowed us to clarify the role of the strong axial and trace anoma- 
lous operators in hadronic weak decays involving c, b or t quark loops. Our results may 
be summarized as follows: 

— The virtual c contribution to K decay amplitudes provides us with an ideal labora- 
tory to study the trace anomaly effects since the scalar gluon density could softly convert 
into a pair of light pions. However, in the Standard Model, the short-distance rise of the 
dimension-seven trace and axial anomalous operators is necessarily 1/M^-suppressed due 
to the chiral structure of the four-quark effective weak interactions, while it is shown that 
there is no trace anomaly effect from dimension-eight operators. Besides, from factor- 
ization arguments, we conjecture a small impact of c quark induced strong anomalies on 
hadronic B decays. Notice that a complementary conclusion has been reached recently by 
the authors of ref.|2l] for the axial 'anomaly tail' due to highly virtual gluons, i.e. above 
the mj, scale. 

— The formalism applies to b quark loop effects in D meson decays too. However, these 
suffer in any case from a multi-Cabibbo suppression. 

— The t quark induced anomalous operators, on the other hand, require to go beyond the 
four-quark approximation. Yet, a simple dimensional argument leads us to conclude that 
the chiral suppression in K{B) decays is now controlled by the m s n^/mt ratio. 

Consequently, if the strong axial and trace anomalies have a role to play in hadronic 
weak decays, in particular in B — ► ij'Ks and K — > tttt, they should arise from either 
nonperturbative effects induced by the light s, d and u quarks or new physics beyond the 
Standard Model. 

Our formalism also provides us with a systematic classification of charm-induced op- 
erators for K decays. In particular, a complete set of gauge-invariant dimension-eight 
operators is displayed. Their effect, of the order of m 2 K /m 2 , might compete with other 
non-leading short-distance QCD corrections |15j . As such, they obviously deserve further 
attention, independently of their anomalous content. Let us notice that distinct l//x 2 
effects beyond dimensional regularization have been discussed in ref . |25j . 

Finally, our method is readily extended to account for c quark contributions to radiative 
K decays. Indeed, the inclusion of photons simply proceeds trough the substitutions: 

9sG^ v — > QsG^v + eqF^ 
9sD ai .. .D an G > QsD ai . ..D an G ^ u + &qd a ^...d OLn F^ LV 

with F^ v , the photon field-strength tensor and q, the heavy quark charge in units of the 
electron charge e. 
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